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Abstract. We analyze operator splitting methods applied to scalar equations 
with a nonlinear advection operator, and a linear (local or nonlocal) diffusion 
operator or a linear dispersion operator. The advection velocity is determined 
from the scalar unknown itself and hence the equations are so-called active 
scalar equations. Examples are provided by the surface quasi-geostrophic and 
aggregation equations. In addition, Burgers-type equations with fractional 
diffusion as well as the KdV and Kawahara equations are covered. Our main 
result is that the Godunov and Strang splitting methods converge with the 
expected rates provided the initial data is sufficiently regular. 

1. Introduction and main results 

We consider operator splitting applied to a class of evolution equations having 
a nonlinear "transport part" , and a linear (local or nonlocal) "diffusion part" or a 
linear "dispersion part". These equations are posed on ~R N for N — 1, 2, 3, and are 
of the form 

u t + div (uv(u)) = A(u), (1.1) 

where u(t, x) is a scalar function and div is the spatial divergence operator. 

The vector-valued operator v( ■ ) and the real- valued operator A{ ■ ) are linear and 



satisfy a set of hypotheses (given in Definition 2.3 below). These hypotheses are 
met, for example, by the popular fractional quasi-geostrophic [3l IH HI QUI [19] 
and aggregation equations [T| 12] UT ] IT2 ] |231 [2^1 125] : 

u t +div(uv(u)) + (-A) a/2 u = 0, v(u) = curl(-A)-' 3/2 w, (1.2) 

u t + div (uv(u)) + {-A) a/2 u = 0, v(u) =\7<f>*u, (1.3) 

where /3,a > 1. Both these equations are paramount examples of so-called active 
scalar equations, that is, equations in which the advection velocity v(u) is deter- 
mined by the scalar unknown u itself. In three dimensions, the general formulation 



(1.1) encompasses also the active scalar equation [5] 

u t + div(uv(u)) — Au = 0, v(u) = divT(u), 

for which well-posedness was established recently in [14] . Here, T is a matrix of 
Calderon-Zygmund operators such that divv(u) = divdivT(u) = 0. This equation 



is a generalized 3D version of the quasi-geostrophic equation (1.2) (a = 2). 

In one dimension, our linearity requirement onn(') limits the type of equations 
we can consider to those of Burgers type, such as 

u t + (u 2 ) x = u xxx , (KdV), 
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Ut + (u 2 ) x = u xx , (viscous Burgers), 

u t + { u2 )x — ~ u xxx + u xxxxxi (Kawahara). 

Since the class of equations studied herein covers a variety of physical models, 
we will give a proper discussion of applications at the end of the paper. 

The main topic of the present paper is analysis of operator splitting methods 



for constructing approximate solutions to (1.1 ). The tag "operator splitting" refers 
to the classical idea of constructing numerical methods for complicated partial 
differential equations by reducing the original equations to a series of equations 
with simpler structure, each of which can be handled by some efficient and tailored 
numerical method. We do not survey the literature on operator splitting here, 
referring the reader instead to the bibliography in |15) . 

The purpose of this paper is to prove, in the context of active scalar equations 



(1.1 1, the well-posedness and convergence rates for two frequently used operator 
splitting methods. Both methods are based on applying repeatedly the transport 
operator and the diffusion/dispersion operator A in separate steps. This splitting 
is very reasonable since one can then use "hyperbolic" numerical methods in the 
transport step and "Fourier space" methods in the diffusion/dispersion step. The 
methods are in the literature referred to as Godunov and Strang splitting and are 
widely used for both numerical computations and analysis. The reader can consult 
[15] for a recent survey of theory and applications; see also for analysis of 
splitting algorithms for the incompressible Navier-Stokes equations. 

Let us now discuss our splitting methods in more detail. For this purpose, we 



first recast (1.1) in the form 

u t = C(u), C{u) = A{u) + B(u), 

where we have introduced the operator B{u) := — div(w(u)). We can then con- 
struct two solution operators <& A and $b associated with the abstract ordinary 
differential equations 

d t <$> A (t, u Q ) = A($ A (t, no)), $4(0, M ) = u 0) 

dt$> B (t,u ) = B($ B (t,u )), $s(0,m o ) = u . 

The first method we will consider, Godunov splitting, is defined as follows. For 
At > given, construct a sequence {u n , u n+1 / 2 }]^i of approximate solutions to 



(1.1) by the following procedure: Let u° — uq and determine inductively 

u n+l ' 2 = $ B (At,u n ), u n+1 = <S> A (At,u n+1 / 2 ), n = 0, ... , [T/At\ - 1, (1.4) 

where \_z\ gives the greatest integer less than or equal to z. 

For Godunov splitting we prove that it is well-posed and that it convergences 
linearly in At. Specifically, we prove the following theorem. 

Theorem 1.1. Let T > be given and assume Uq € H k with 6 < k e N. Then, 
for At > sufficiently small we have the following: 



(1) The Godunov method (1.4) is well-defined with 

IKH^<C, n=l,...,[T/At\. 

(2) The error satisfies 

|K - u(nAt)\\ Hk -^ x{a , 2i < C||tio|| H »At, (1.5) 

where a is the highest number of derivatives occuring in A. 

The other method we consider, Strang splitting, is defined as follows: For At > 
given, construct a sequence {u™, w n+1 / 4 , u™ +3 / 4 }| l ^j' 1 A ^ of approximate solutions 
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to ( fTTTj ) by the following procedure: Let u° = uo and determine inductively, for 
n = 0, . . . , [T/At\ - 1, 

U «+V4 = $ B (i At, u"), u"+ 3 / 4 = $.4 (At, u" +1 / 4 ), 

? (1.6) 
M n+1 = $ s (-At,u™ +3/4 ). 

For the Strang splitting algorithm we prove well-posedness and second order 
convergence, provided the initial data are sufficiently regular. 

Theorem 1.2. Let T > be given and assume Uq € H k with 6 < k e N. Then, 
for At > sufficiently small we have the following: 



(1) The Strang method (1.6) is well-defined with 

\\u n \\ H ><C, n=l,...,[T/At\. 

(2) The error satisfies ( |1.5[ ) and 

\\u n - u(nAt)\\ Hk - 3 ^ {aA} < C{At) 2 . (1.7) 



The a occurring in (1.7) is the highest number of derivatives in A. 



The approach leading to Theorems |1.1| and |1.2| utilizes the analysis framework 
put forth in the recent paper |16j for the KdV equation. This framework works 
with a definition of the Godunov method given in terms of a specific extension of 
the splitting solution {u n , u n+1 / 2 } n to all of [0, T]. An adaption of this extension 
that also applies to Strang splitting was provided in 17J. This extension, which 
is different from the one used in [TH] for Strang splitting, will be employed herein. 
This allows us to treat in a unified manner a rather general class of equations, 
covering those treated in [IB] , [IT] , and [TS] , as well many additional equations not 
treated in these papers. At variance with [17], we do not require the divergence 
of the velocity field v(-) to be zero, thereby enlarging significantly the class of 
equations that can be handled. 

The remaining part of this paper is organized as follows: Section [2] collects some 
preliminary results needed later on. The Godunov and Strang splitting methods 
are analyzed in Sections [3] and [4j respectively. 

2. Preliminary existence and regularity results 

2.1. Notation and two technical estimates. We will use LP to denote the 
Lebesgue space of integrable function on M. N with exponent p. If I denotes a N- 
dimensional multi-index, i.e., I = (li, . . . , In), lj € No, we write 

to denote any derivative of the Ith. order. If I £ N, we let 

V'f = {V l f\\l\=l} 

and 

VV:V^=^VVV' 5 . 

\i\<e 

We will be working the Sobolev spaces 

H k = {fes' | (i + \d\ 2 ) k/2 Hf(0)eL 2 }, 
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(where S' denotes the set of tempered distributions) and T denotes Fourier trans- 
form. If A; is a natural number, H k is the standard Sobolev space with inner product 
and norm given by 

k 

(f,9) H > =E<vV,V^)i», \\f\\ H * = (fJ))£, 



e=o 



where we have introduced 



\l\=t 

Throughout the paper we will strongly rely on the following two technical lem- 
mas. The validity of these estimates are the primary reasons for the requirements on 
v given in Definition 2.3 Their proofs are straightforward, but somewhat tedious. 
For this reason, proofs are deferred to the appendix. 

Lemma 2.1. Let k > 6. Then 
k 



E 



V s (div(/u(/))) : V s / dx 



< C\\f\\ Hh -4f\\ 2 Hh , feH k 



Lemma 2.2. Let k > 4. Then the following estimates hold 
k 

1 2 f „ rrk 



£ / V s div(fv(g)):V s f dx 
s =o Jrn 

V [ V s div(gv(f)):V s f dx 

tit JS- N 



<C\\g\\ H -\\f\\H*, f,9&H K , 



<C\\g\\ Hk+ i\\ff Hk , feH k ,geH 



k+l 



(2.1) 



(2.2) 



2.2. Existence and regularity results. Our splitting methods are based on al- 
ternately solving the equations 

U t =A(u), u\ t=0 = u , (2.3) 

u t + div (uv(u)) = 0, u\ t=Q = u - (2.4) 
To analyze the methods we will need to have some control on the behavior of 
solutions to each of the two equations separately, at least on short time intervals. 
Before we start discussing available results, let us state the specific properties we 
will require of the operators v and A. 

Definition 2.3. We say that the operators v, A are admissible provided: 

(1) A and v are linear operators satisfying the commutative property 

A(vi{'))=Vi(A(-)), i = l,...,N, 

where is the ith component of v. 

(2) v : L p — > L p is bounded for any p < oo, and, if N > 2, 

||divt>(u)||LP<C||«||L„ N>2. 

(3) A is a differential operator; specifically, there is a positive integer a such 
that A: W a+k ' p — > W k,p , is bounded for all k < oo and p < oo. 

(4) A is either conservative or diffusive: 

A(u)u dx<0, ue W a+k > p . 

(5) A satisfies the following commutator estimate 
\\A(fg)-fA(g)-gA(f)\\ Hk < C\\f\\ H u^^y-i IMI^w^w, 
for all f,ge H k+ma ^ a ^- 1 for k > 3. 
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(6) The equation ( 1.1 1 is well-posed in the sense that for any given finite T > 
and initial data uo G H k , there is a solution of satisfying 

u G C([0,T];H k ). 



The first equation (2.3) is a linear equation with constant coefficients ((1) in 
Definition 2.3) that preserves (or diffuses) all Sobolcv norms over time. From this, 
the following lemma follows readily. 

Lemma 2.4. // uo G H k for some a < k < oo, then there is a unique solution 
ueC([0,T};H k ) of (pi. 



The nonlinear equation (2.4) is of hyperbolic type. Typically, the best available 
existence results for this type of equations are local in time existence of smooth 
solutions. 

Lemma 2.5. Let uq G H k for some k > 4. There exists a time T > and a 
function u G C([0,T); H k ) D C 1 ([0, T); H k ~ 1 ) such that u is a unique solution of 

(2.4) on [0,T). Moreover, if T is the maximal time of existence for (2.4), then 

\\m\\u{t)\\ Hk = oo. 

Proof. We will prove the local in time existence by demonstrating compactness of 
solutions to the approximation scheme: 

uf + b^ 1 ) • Vu m + u m div Bfu™- 1 ) = 0, w m | t=0 = u . (2.5) 

For this purpose, we let u° = u and sequentially determine the sequence {u m }^ =1 
as the solutions to the app roxi mation scheme (2.5). Now, for each given tt m_1 , the 
requirements of Definition [2~3] yields ^(w™" 1 ), divu^" 1 " 1 ) G C([0, T];H k ). Thus, 

( 2.5 ) is a linear transport equation with smooth coefficients and consequently admits 
a smooth solution for all times. 

Let us now calculate the H k norm of u m . To achieve this, we apply V s to (2.5), 
multiply with V s u m , sum over all s = 0, . . . , k, and integrate to obtain 

dth\u m (t)f Hk = - V / (V s div(u m v(u m ^)) : V s u m )(x,t) dx 

<C\\u m -\t)\\ Hk \\u m (t)f Hk , 

where the last inequality is an application of Lemma |2.2| Applying the Gronwall 
inequality to the previous inequality yields 



Next, we fix 



IK 1 Hl~([cw) < e CT||, ' m " 1|l ^([P.n i «*)||« || Hfc . 
log 2 



(2.6) 



T < 



2C\\u \\ Hk 



where the constant C is the one appearing in (2.6). Using the bound on T in (2.6) 
and iterating the resulting inequality (starting from m — 1), we obtain 

ll U " l ||L~([0,T];ff fc ) < 2||uo||fffc- 

Next, let us calculate the H k ~ 1 norm of it™. By direct calculation, 
k-i . 

\\u?{t)\\ H u-r = / (V s div(u m ^( M m - 1 )) : V s u?){x,t) dx 

< Cr||« TO (*)]| jH - fc ]|« TO - 1 (*)|U* 
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At this point, we have proved that u m E C([0, T); H k ) n C 1 ([0, T); H k ~ x ), inde- 
pendently of to. Hence, we can assert the existence of a function u <E C([0, T); H k )D 
C 1 ([0,T);H k - 1 ) such that 

u m u in C([0, T); H fc ) n C x ([0, T); iJ^ 1 ), 

as to — > oo, where the convergence might take place along a subsequence. Compact 
Sobolev embedding then tells us that 

u m ^u, mC([0,T);H k - 1 ), 

again along a subsequence. Note that this is not sufficient to pass to the limit in 
(2.5). Indeed, we need that the whole sequence converges. To prove this, we let 
w m = u m — a" 1 " 1 and observe that 

w™ + div(^(w m - 1 )w; m ) + div(«(u; m - 1 )u m ~ 1 ) = 0. 



It follows that 



-t4|K1^ =- [ l\w m \ 2 divviu" 1 - 1 ) + div (viw™- 1 ^" 1 - 1 ) w m dx 



< C\\u m - x \\ Hk (|| 



w m \\ 2 L2 + \\w™- l \\l 2 ) 



where we have used the requirements on v(-) (Definition 2.3), Sobolev embedding, 
and the Cauchy inequality. An application of the Gronwall inequality, using that 
w m (0) = 0, we obtain 



\\w m \\Mt) < Ce CT j Wur-^h' ds < ± ^- sup l^ 1 



(tCe CT Y 

j »Lip \\W \\ L 2 

m - te(o,T) 



rim 

< j" SUP || W || £,2 — > 0. 



From this we conclude that u m is a Cauchy sequence in L 2 and hence that the entire 
sequence converges. We can now pass to the limit in (2.5 1 to conclude the existence 
part of the lemma. Uniqueness is an immediate consequence of the regularity of 
the solution (i.e subtracting one solution from a possibly different solution, adding 
and subtracting, and using the regularity). 

The blow-up at maximal time of existence follows from u being continuous in 
time. □ 



3. GODUNOV SPLITTING (PROOF OF THEOREM 1 1 - 1 [ ) 

To prove Theorem we will utilize the analysis framework put forth in |16j . 
To this end, we introduce an extension of the splitting solution {u 11 , u n+1 / 2 } n to all 
of [0, T]. The definition is posed on the two-dimensional time domain 

[T/AtJ-l 

^At = [J [tn,t n +l] X [tn,tn+l], 

and goes as follows: 

Definition 3.1 (Godunov splitting). For At > given, we say that -d is the 
Godunov splitting approximation to (1.1) provided that $ satisfies 

0(0,0) = O , 
0«(*>*n) = B(ti(t,t n )), te(t n ,t n+ i], 

T (*> t) = A(0(t, r)), (t, t) € [t n , t n+1 ] x {t n ,t n+1 }. 
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Observe that 

tf(t n ,t n )=u n , n = 0, . . . , [T/At\ -1. 

Thus, i?(t, t) is indeed an extension of u n to all of [0, T]. 
To measure the error, we will use the function 

e(t) = t) - u(t), 

where u is the (smooth) solution of 

Since i? is an extension of u n to all of [0, T], it is clear that Theorem is an 
immediate consequence of the following lemma, which is our main result in this 
subsection. 

Lemma 3.2. Let T > be given and assume Uq £ H k with 6 < k € N. Then, for 
At sufficiently small we have the following: 



(1) The Godunov method in Definition 3.1 is well-posed with i? € C(£l^ t ]H k ). 

(2) The error e(t) satisfies 

\\e{t)\\ Hk - m » xia , 2} = \\$(t,t)~u(t)\\ Hk - n ^ {a , 2i < C\\u \\ 2 Hk tAt. 



The proof of Lemma |3.2| will be a consequence of the results stated and proved 
in the ensuing subsections. The closing arguments will be given in Section [3. 3| 



3.1. Error evolution equations. The main benefit of the new definition of the 
Godunov method is that it allows us to derive continuous-in-time evolution equa- 
tions for the error e. Though these equations was previously derived in [TB], we 
repeat the derivation here for the convenience of the reader. For this purpose, we 
need the following Taylor expansion that holds for any smooth operator E 

E(f + g) = E(f) + dE(f)[g] + f\l - j)d 2 E(f + lg )[gf d r 

Jo 

Using the definition of $ and the above Taylor formula, we deduce 
e t - dC(u)[e] = <& t + <& T - u t - dA(u)[e] - dB(u)[e] 

= d t + A(<&) -(A + B){u) - dA{u)[e] - dB(u)[e] 
= 6 t - B(6) + - A(u) - dA(u)[e}) 
+ (B(&) - B(u) - dB{u)[e\) 

= F(t,t)+ f\l-j)d 2 C(u + 1 e)[e} 2 d 7 , 



(3.1) 



where we have introduced the "forcing" term 

F(t,T) = Mt,r)-B(#(t,T)). 

By direct calculation, 

F T - dA($)[F] = v tT - B{d) T - «L4(0)[0 t ~ B(d)] 

= A(#) t - dB($)[tf T ] - dA(tf)[0 t ] + dA{tf)[B(tf)} 
= dA(d)[i)t} - dB(6)[A(0)] - dA{$)[$ t ] + dA(ti)[B(ti)] 
= [A,B}(#), 

where we have defined the commutator 

[A,B](f)=dA(f)[B(f)]-dB(f)[A(f)]. 

In our case, the operator A is linear and hence 

dA(f)[g}=A(g), d 2 A(f)[g,h}=0, (3.3) 



(3.2) 
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while the operator B satisfies 

B(/) = -div(.M/)), 
dB{f)[g] = - dW(fv(g) + gv(f)), (3.4) 
d 2 B(f)[g,h] = -dbr(hv(g)+gv(h)). 
This and the requirement (1) in Definition |2.3| yields 

[A, B](f) = -A (div(/«(/))) + div (fv(A(f)) + A(f)v(f)) 
= - div (A(fv(f)) - fA(v(f)) - A(f)v(f)) , 

which is exactly the divergence of the commutator appearing in (5) of Definition 
|2.3| Thus, the following lemma follows directly from this requirement. 



Lemma 3.3. Let f £ H , where a is given by (3) in Definition 2.3 Then, 
\\[A,B](f)\\ Hk ^ {2 , a} <C\\f\\ 2 Hk - 



Using (3.3) and (3.4) in (3.1) and (3.2), we obtain the following evolution equa- 
tions for the error e 

e t + div (ev(e) + uv(e) + ev{u)) - A(e) = F, t e (0,T) (3.5) 

F T - A(F) — [A, (t,r)Gn At . (3.6) 

Upon inspection of these equations, we see that the error e satisfies an equation 
similar to (1.1), but with an additional source term F. 

3.2. Estimates valid under the assumption of regularity. In this subsection 
we state and prove some results that will be needed in order to prove Lemma |3.2| 
Due to the special configuration of the time domain, it will be convenient to use 
the following notation for all times prior to a given time (a, Q) £ i^At : 

n a ^ t = {(t, t) £ f2 A t I o < t < a, < r < c} , 
The following lemma is the most essential ingredient in the proof of Lemma |3.2| 



< 



Lemma 3.4. Let i9(0,0) := Uq £ H k . Assume the existence of a time (a, £) £ 
[0,T] 2 and a finite constant 7 > such that 

(i,r)£^t C - 



\W(t,T)\\ Hh -» 



t { Q ,2} 



There is a constant C(~f) determined by 7, uq, and T such that 

||0(t,T)|U*<C(7), (t,T)£0 A i 

Proof. Fix any (t,r) £ ^ A J" and let n be such that t,T £ [t n ,t n+ {\. Using the 
definition of d (Definition 3.1 1, we see that 

i 1 
dr-\\d{t,T)\\ 2 He = V / V S ?9 T : V s i? dx = V / A(V S ^) : \7 S $ dx < 0, 
2 ' S=0 A« s= o«" 



for any I = 1, . . . , k, where the last inequality is requirement (4) in Definition 2.3 
The previous inequality yields 



¥(t,T)\\ Hk < \\0{t,t n )\\ H u. 



(3.7) 



Next, we let i > 4 be an integer and apply Definition 3.1 to obtain 

1 

dthm,t n )\\ 2 He = -V f V s div(#v(#))--V s #dx 

1 ~^JS- N 



<c\\n H i-4m t , 
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where we have applied Lemma 2.1 If £ is such that £ — 2 < k — max{2, a}, we can 
apply the Gronwall lemma to the previous inequality to obtain 



m,t n )\\ H * < \m n ,t n )\\ He e°^ At . 



(3.8) 



By combining ( |3.7[ ) and (3.8), we thus obtain the bound 

\\&(t,T)\\ H i < ||u ||^e C7T , 4< i< fc-max{2,a}. (3.9) 

Now, we can repeat the above arguments with a new 7 := C(-~f) = \\uo\\jjte tT , 
to obtain 

||i?(t,r)||^ +2 < \\uo\\ H t+ a e a ^ T , 4 < £ < k — max{2, a}. 

Clearly, we can repeat this process n times until £ + In = k (a even) or until 
£ + 2n = k — 1 (a odd). In the last case, we can conclude that 



\\i9(t, T )\\ Hk - 2 < H^t)!!^-! < || Mo ||^-ie c ^ T , 



and hence the arguments leading to (3.9) hold for / = k. 



a 



In the following lemma, we give our main error estimate. Note that the result 
is only valid under a strong regularity assumption on the splitting solution 1?. This 
assumption will be fully justified when we prove Lemma 3.2 in the next subsection. 

Lemma 3.5. Assume the existence of (a, Q € [0, T] 2 such that the splitting solution 



Then, 



||0(t,T)|| JJ *<<7( 7 ) 1 (t,r)en2i. 



|e(t)|| Hfe - m « { «,2 } <C(7)At, t<a. 



Proof. To shorten the notation in this proof we write £ = k — max{a, 2}. Let 
us commence by estimating the size of the source term F. For this purpose, we 
apply V s to (3.6), multiply the result with V S F, integrate by parts, and sum over 
s = 0, . . . , k, to obtain 

1 ' r 

-d T \\F(t, T )\\ 2 He =Y / (A{V S F) : V S F + V S {A,B}($) : V s F) dx 

2 ' I^k~ 



< 



V / \7 s {A,B](tf) : V S F dx 



<C\\F\\ Hl \\[A,B]W\\ Ht 



Since £ = k — max{a, 2}, we can apply Lemma |3.3| to obtain 

d T \\F(t,r)\\ H i < \\[A t B](#)\\ Hi < C\\d\\ 2 Hk < C( 7 ). 

By definition F(t, t n ) =0. Hence, integration in the t direction from t n to r yields 

\\F(t,T)\\ H i<C(i)M. (3.10) 

Let us now turn to the evolution equation ( |3.5[ ) for the error. By applying V s 
to (3.5), multiplying the result with V s e, integrating by parts, and summing over 
s = 0, . . . , k, we obtain 

ld t \\e(t)f Hl =Y) f / A(V s e) : V s e dx + [ V s F:V s edx) 

t^o u ' Rjv "' RJV 

- 2~J / VS ( div ( ev ( e ) + uv ( e ) + ev i u ))) '■ dx 

s=0 jRN 

< C\\F(t,t)\\ Hl \\e(t)\\ Ht +C\\u(t)\\ HW \\e(t)\\ 2 Ht , 



10 



HOLDEN, KARLSEN, AND KARPER 



where the last inequality is an application of Lemma 2.2 and Holder's inequality. 
Now, since uq £ H k , requirement (6) in Definition 2.3 tells us that the analytical 



solution u £ C([0, T\; H k ). Thus, the previous inequality leads us to the conclusion 

dt\\e(t)\\ H t < C\\F(t,t)\\ H i + C\\e(t)\\ B , < C (At + \\e(t)\\ B i) , 

where the last inequality is | |3.10 1. An application of the Gronwall inequality con- 
cludes the proof. □ 



3.3. Proof of Lemma |3.2[ We will determine the maximal size of At > during 
the course of the proof. Specifically, the size of At will be determined in accordance 
to the maximal existence time of solutions to the nonlinear transport equation ( 2.4 1. 
For the convenience of the reader we recall that the maximal existence time T* of 
solutions u to (2.4) is characterized by 

(3.11) 



lim 



\U\\ H k = 



1. The proof of well-posedness ((1) in Lemma 3.2) will be a direct consequence 
of the following lemma. 

Lemma 3.6. Let k > 6 and uq £ H k . There exists constants ft, 7, depending only 
on uq and T , such that if At < ft and 

\\d[t,T)\\ Hk -^ {a ^ < 7 , (t,r) £ Q^i 
for some (a, £) £ flAt, then 

7 u _^ n *X 



\W,T)\\ H >- 



< 



(t, r) £ n 



At ■ 



Let us for the moment take Lemma 3.6 for granted and explain why it concludes 
our proof of L emma |3.2| For this purpose, we let At < ft, where ft is dictated 
by Lemma 3.6 Since the initial data $(0,0) = u £ H k , our short-time existence 



result (Lemma 2.5) enables us to start constructing $(£, 0) according to Definition 



3.1 up to some time T* , that is, to start applying the B operator (the nonlinear 



transport equation (2.4)) in the first time step. Now, on the interval (0, T*) there 
is no problem determining 7 in accordance to both Lemma |3.6| and such that 

\\ti{t,0)\\ Hk -^ x{a ,2 } < 7 , 0<t<T*. 



Lemma 3.6 can then be applied to conclude that 



\m,o)\\H>- 



< 



7 



< t < T* 



However, in view of (3.11), and since the norm is continuous in time there must 
exist an e > such that 

mt,0)\\ Hk - mia , 2} < 7 , 0<t<T*+e. 

Thus, by repeating these three steps, we can conclude that 

||tf(i,0)|| fffc - max{Q , 2} < 7 , 0<t<At. 



We can now start to determine i9(t,r), for < t < At, according to Definition 3.1 



that is, to start solving with the A operator (the linear equation (2.3)). By virtue 
of our existence result for (2.3) (Lemma 2.4), we can repeat the above argument to 
conclude that $(t, t) satisfies 

\\<&{t, T)|| fffc - m ax {Q , 2} < 7, < t, t < At. 

Besides well-posedness of the first time step, this result tell us that 7 also bounds 
the # fc - max { Q > 2 } norm of i?(At, At), which is indeed the initial data in the next 
time step. Hence, we can repeat the entire process for all time steps recursively to 
conclude that 

\\0(t,T)\\ H >- mi *n <7, (t,r)en A t. 
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Lemma 3.4 can the be applied to obtain 

\m,r)\\ Hk <C(j), (t,r)en A u 



which concludes the proof of well-posedness, i.e., (1) of Lemma 3.2 



2. Since we now know that ||i?(t,T)||jj* < C(7), (t, r) £ &>At, the postulates of 
Lemma [3~5| are satisfied. Thus, 



K*)lliT*- 



,. 2} <AiC, i€(0,T), 



which is precisely (2) in Lemma 3.2 



We conclude this section by proving Lemma |3. 6 



Proof of Lemma \3.(\ To shorten the notation, we let t — k — max{a, 2}. Assume 
that (a, C) £ SIa* is such that 



||i?(t,r)|| H i<7, (t,^^, 

where the values of 7 and At are still to be determined. 
From Lemma 3.4 we know that ( 3.12[ ) implies that 

\m,r)\\ Hk <C( 7 ), (t,r) efi& 

We can also apply Lemma [3~5| to obtain 

\\e(t)\\ H e <AtC( 7 ). 



(3.12) 



(3.13) 



Now, fix any (t, r) £ . By definition of the Godunov method (Definition 3.1 ), 

k k 



dr^\W(t,r)\\ 2 H , 



V / V s i? r : V s i? dx = V f A(S7 S $) : V s tf dx 



< C\\d\\ B t\\A(&)\\ H i < C\\4\\ H 44\\ Hl+a < C\\0\\ H 4#\\ H >, 
from which it follows that 

\d T \\d(t,T)\\ Ht \ <c\\n Hk . 

Using this inequality and Lemma |3.4[ we find 



\#{t,T)\\ H t < \\4(t,t)\\ H t + J \d T \\ti(t,f)\\ Hk \ d\Ulder 

< \\0(t,t)\\ Ht + AtC sup ||0(t,f)|| H » 

fe[t,T] 

< \\0(t,t)\\ H *+AtC(i) 

< ||e(t)|| ff * + ||u(t)||^+AtC(7) 

< AtC(-f) + C 2 , 



by applying (3.13), we obtain the estimate 

||i?(t,r)|| ff i < ||e(*)|U/ + ||«(t)|| ff i+AtC( 7 ) 
< AtC{ 1 ) + C 2 . 



(3.14) 



Finally, by fixing 7 and At according to 



7 - 4C 2 , 



C, 



At < -^ 



C(7) 



\W,t)\\h* < 



7 



and applying this information to (3.14), we see that 
which concludes the proof. 



□ 
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4. Strang splitting (proof of Theorem 1.2) 

As for the Godunov method, our convergence analysis will require a continuous 
definition of the Strang method. In contrast to the Godunov case, we will now 
introduce three time variables instead of two [T7] . We consider the domain: 



[T/AtJ-1 

U i 

n=0 



y. — ] x [*n»*n+l] x Ly — ] 



The continuous Strang method is given by the following definition. 

Definition 4.1. For At > given, we say that is the Strang splitting approxi- 
mation to (1.1) whenever solves 



0(0,0,0) =«o, 

**(*>*», |) = B(0(*,t»,|)), *e(|,^], 

T (i,r,^) = A(0(t,r,|)), (t,r) G x (t n ,t n+1 \, 

K(t, r, u) = B(fi(t, r, U )), (t, r, u) G x t„ +1 ] x (|, 



(4.1) 



In each box f^*- , 



i+iJ x ~ 2 2 li ]i we wm mainly consider the function 
along the diagonal, i.e., the function 0(|,f, |) for t G [t n ,t n+ i]. Observe that 
each point on this diagonal is a Strang splitting solution for a specific time step. 
More precisely, 0(f,t, §) is a Strang splitting solution with time step t ~ t n . An 
easy consequence is that 



n=0,...,LT/AtJ, 



and hence that 0(|, £, |) can be seen as an extension of {u n } n to all of [0,T]. 
To measure the error, we will use the function 



e(*)=^»*»|) -«(*). 



where u is the (smooth) solution of 

Theorem |1.2| is an immediate consequence of the two following lemmas. 



Lemma 4.2 ( Well -pos edness) . Let the Strang splitting solution of (1.1) in the 
sense of Definition 4.1 and (4.1). Let T > be given and assume uq £ H k with 
6 < k G N. Then, for At > sufficiently small 

(1) TTie Strang method is well-posed with G C(£lAt', H k ). 

(2) TTie error e(i) satisfies 

||e(t)|| ff *-„« { ..»} = ||i?(|,t,^)-u(t)||B*-»«{-.»} <*CAi. 



Lemma 4.3 (Con verg ence) . Let f/ie Strang splitting solution of (1.1) m i/ie 
sense of Definition 4-1 and (4.1). Let T > be given and assume uo G -ff fe with 
6 < k £ N. Then, for At > sufficiently small 

\\e(t)\\ Hk - Sm ^ aA} < C(At) 2 . 



Lemmas 4.2 and 4.3 will be consequences of the results stated and proved in the 
ensuing subsections. 
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4.1. Error evolution equations. By direct calculation, we see that the error e 
satisfies the time-evolution 

e t - dC(u)[e] = y + ?? T + y - u t - dA(u)[e] - dB{u)[e] 

= y +# T + \b(Q) -(A + B){u) - dA{u)[e] - dB{u)[e] 

+ (A(tf) - A(u) - dA(u)[e\) + - B{u) - dB(u)[e]) 

= F(t)+ [ (l-7)d 2 C(u + 7 e)[e] 2 dry, 
Jo 

where F(t) = F(f,t, |) and 

F(t, r, w) = i (tf t (t, r, w) - B(#(t, t, w))) + ?9 r - A(tf (t, r, w)). 

Since ^ r — A(#(i, r, w)) = when ui = if, n = 0, . . . , [F/AtJ — 1, we can apply 
the arguments of ( |3.2[ ) to obtain 

^ r - dA(0) [F] = i [A, B] (0), (t, a, y ) G n At 

where as before 

[A, J B](/)=dA(/)[B(/)]-dB(/)[A(/)]. 
We also derive the following equation for the evolution of f in w: 
F u -dB(#)[F] 



dA(<ff)[0 u ]-dB(<ff)[9 T -A(0)] 

= \B{*) t ~ \dB(W»\ - \dB{mt] + \dB{d)[B{d)] 
+ dB(#)[0 r ] - dA($)[B(-d)} - dB(i9)[$ T ] + dB(&)[A(d)] 

= i (dB(0)[0 t ] - dB(0)[B(0)] - dB($){$ t ] + dB(0)[B(0)]) 
+ dB(tf)[A(tf)] - dA(0)[fl(i?)]. 
Thus, recalling the definition of [ • , • ], we find 

F u -dB(fi)[F] = [B,A](ti). 



For the class of equations ( 1.1 1, the evolution equations read: 

e t + div (ev(e) + uv(e) + ev(u)) — A(e) = F, t e (0,T) (4.2) 

F T -A{F)= l -[A,B]{d), (t,T,^)en At , (4.3) 

F w +div(lMF) + Ft7(0)) = -[A,B](0), (t,T,w)efiAt. (4-4) 

As for the Godunov method, it will be convenient to define a notation for all 
times prior to a given time (a, (, v) G flAt- 

= {(t,T,u)en At \o<t<a, o < t < c, o < u < v} . 

To prove Lemma |4.2[ we will need Strang versions of Lemmas |3.4| and |3.5| 



Clearly, there is no problem to extend Lemma 3.4 to conclude the following re- 
sult. 
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Lemma 4.4. Let ??(0,0) := uq G H k . Assume the existence of a time (<r, r, v) G 
[0, T] 3 and a finite constant 7 > such that 

\\ti(t,T,u)\\ H k-^ { = ,2} < 7, (t,T,Uj) G ^At C,1/ - 

T/ien i/iere is a constant C(j) determined by 7, uq 7 and T such that 



||i?(* > r,w)|| H *<C(7), (t,T,u)6ff 



At ' ■ 



Let us now prove the Strang version of Lemma 3.5 



Lemma 4.5. Assume the existence of (a, £, t") € J7a* suc/i t/iai the splitting solution 



II^t^jiIh* <c(rf), (t,T,u)£n° At 



Ther 



\\e(t)\\ Hk -^ { o,,2 } < A*C(7) > t< er. 

Proof. Let £ — k ~ max{a, 2}. As in the proof of Lemma |3~5j we first estimate the 
size of the source term F. For this purpose, fix any (t, r, oS) G an< ^ ^ n ^ e 

such that (i, r, w) G if, ^r 1 x [t n ,t n+1 ] X t -f, t -^ L . 

Now, in the plane given by u = we have that F satisfies the equation (4.3 1. 
Since this equation is similar to (3.6), we can repeat the arguments in the proof of 
Lemma [3~5l to conclude that 



|F(t,a,|)|| ff ,<AtC(7). 



(4.5) 



Hence, to estimate F at the point (t,T,ui), we can integrate from (t, t, lj) in the to 
direction to the plane given by to = and apply (4.5). To achieve this, we first 
apply V s to (4.4), multiply with V S F, integrate, and sum over s — 0,...,k, to 
obtain 



d u ±\\F{t,T,u>) 



2 



V / V s div (§v(F) + Fv{6)) : V S F dx 

D 

[ V S [A,B}{$) :V S F dx 

J-R N 



E 

s=0 



3.3 Since d G -ff fc , we can 



In the last inequality we have applied Lemmas 
integrate the last inequality to obtain 



2.2 



and 



\\F(t, T ,u)\\ Ht < \\F(t,v,f)\\ H i+C(>y)At<AtC(>Y), 



where the last inequality is (4.5). 



Since (4.2) is identical to (3.5), we can repeat the arguments in the proof of 
Lemma |3.5| to conclude the proof. □ 

4.2. Proof of well-posedness (Lemma |4.2[ ). We can apply similar arguments 
as those of Section 3.3 (for the Godnov method) to prove Lemma 4.2 In particular, 

is a 



upon inspection of the arguments in Section |3.3| it is clear that Lemma 4.2 
consequence of the following result. 

Lemma 4.6. Let k > 6 and i9(0,0, 0) := uq € H k . There exists constants j3 and 
7, depending only on T and u , such that if At < /3 and 



\\&(t,T, Cj)\\ Hk - 

for some (a, £, v) G Oaj, then 



< {Q , 2} <7, (t,r,w)en? c '" 



'At 



\\m.T.J\\\ w . ,„. < J, (t 1 T, U )6^. 



OPERATOR SPLITTING FOR ACTIVE SCALAR EQUATIONS 



15 



Proof. Let £ = k — max{a, 2} and assume that (t, r, oj) € f2^£ ' v is such that 



\\d(t,T,Lj)\\ H , < 7 , (t,T,u)6fl 



At > 



where the values of 7 and At are to be determined. Lemma 4.4 can then be applied 
and yields 

\\ti(t,T,uj)\\ Hk <c( 7 ), v(t,T,u)en2i'°. 

Furthermore, we can apply Lemma |4.5| to conclude the estimate 



||e(t)|| ff *<AtC( 7 ). 

Now, to proceed we fix any (t, r, oj) £ ^At • Using the definition of the Strang 
method (Definition 4.1 1 and Lemma 2.1 we see that 



dJ-\\d(t,r,oj)\\ 2 He 



E 



V s *, : V s i? dx 



1 

V / V s div(tft>(tf)) : V s i? dx 

<c\\0\\ Ht -40\\% t <c\\of He <c 7 3 . 



From this it easily follows that 

\Mt,T,e>)\\ Ht < \\{}(t,r,t)\\ He + Atcij). 

By the same calculations as those in the proof of Lemma |3.6[ we also deduce 
\m,T,t)\\ H * < ||l?(i,2t,t)|| H< +A*C(7) 

< ||e(2t)|| ff « + ||u(2t)|| ff «+AtC(7) 

< C{j)At + C 2 . 



(4.6) 



Using this in (4.6) allow us to conclude 

\\0(t,T,w)\\ Ht <C(j)At + C 2 . 



(4.7) 



Since (t, t, oj) was arbitrary, we can conclude that (4.7) holds for all (t, r, oj) € fi^lt 
Finally, we fix 7 and At according to 

C 2 



7 = 4C 2 



At < 



C(7) 



in (4.7) to obtain 



a 



Equipped with the previous lemma, Lemma 4.2 can be proved as we did with 
Lemma 13.21 in the Godunov case. 



4.3. Temporal regularity. To prove second-order convergence (Lemma 4.3), we 
will need regularity in time of our splitting solution. Since we are working with 
three time variables, it will be convenient to recall the notation 



|/| = h + h + h, 



dt^dT^dui 1 ^ 

for a multi-index I = (Z 1; l 2l £3). We will also use the notation 

0111/ 



V t K / = 



dt^dr^du 1 * 



= k 



for any natural number k 
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Lemma 4.7. Let ■& be the Strang splitting solution in the sense of Definition 4-1 
Ifu € H k , then 



|V£tf(i,T,w)|| fffe - mm a* { „,i} <C, (t,r,u) e n 



At- 



Proof. We will argue by induction on \l\. For |Z| = 0, the result follows from Lemma 



4.2 To proceed we assume that the result holds for |/| = 0, ... ,q. To close the 



induction argument it remains to prove the result for \l\ = q + 1. 
To simplify notation, we let 

X = k — \l\ max{a, 1} = A; — (q + 1) max{a, 1}. 

Now, let (t',T',oj') £ f^At be arbitrary and fix n such that (t',T',Ld') e [ 



t„ tn+11 



+ U 



L 2 ' 2 



An arbitrary component of Vjz9 can be written in the form 



To estimate the arbitrary component at the point (i',r',a/), we apply the funda- 
mental theorem of calculus to obtain 



E 



S=0 "2 



V 5 9 w e^(i',T',a;) : V s e^(t',r',a;) dxdCo 
2 ii - 9 



5iie5i(fy,^)ii^ 



A i 



EEEE 



n / \m 



s=0 r=0 n=0m=0 



X 



(4.8) 



where we used the definition of and i? w (cf. (4.1 )) to conclude the last equality. 
Let us consider three separate cases of to + n + r in the quadruple sum above, 
(i) If to + n + r = q + 1, the corresponding term in the above reads 



E 



S=0 "2 



V s div(e«+ 1 £ (t',r',w)- i ;(^',r',^))) : VBg^/^) cfedu) 



<c/ K^fy,*))^ eg-^fy.c) 



da; 



< C 



q8J(*v,*) 



da), 



(4.9) 



where we have applied Lemmas |2.2| and |4.4| 

(ii) If to + n + r = 0, we can also apply Lemmas 2.2 and 4.4 to conclude 



s=0 17 ~ 



V s di 
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< c I \\wy,a)\\ H x + 4Qj£tfy,Q)\\ 3 H > <& 

/oj' 



,ui)\\ 2 h x du. 



(iii) For the remaining cases (1 < m + n + r < q), we apply the Holder inequality 
to obtain 



Er / v-div(e^(^y,c)«(e«+^:-(t'y,a;); 

<C T 1 1 (f , ^' , c^) 1 1 1 1 ©^+^" (i' , ^' , cD) 1 1 (4.10) 

J if- 

x||e^T:--:™(*',r',a})|| H x +1 ^ 



where we have used that our induction hypothesis yields 

||V?i?|Ua +1 = ||V?i?|| ff » + i-(, + i)»« { a.i> < HV^Hff*-,^-,!} < c. 



By applying (4.9)— (4.10) to (4.8), we gather 



\\m\{t>y,j)\ 



2" ^ 



2 



< 

- 2 



+ c iie^^y.^ii^ds 



where the last inequality as an application of the Holder inequality to the second 
term. Now, by applying the Gronwall inequality to the previous inequality, we 
conclude 



(4.11) 



We have now derived a bound on 9?~t^(t',T',u/) in terms of Q^g(t', r', ^f). 
Next, we derive a bound on Q?t^(* / ,T / , %-) in terms of 8^^(t',i„, %*>). For this 
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purpose, we once more apply the fundamental theorem of calculus to obtain 

£lieSJ(<v.y)llir* 



1,3 
1 

2 11-/,; 
A 



t r, 



' s, -^) dxds 



2 ll®i,j,*(*'> "oOllff* 



<5fl©SJ('.*»>|)lll*. 

where we have used that d T — A(i3) for uj — and (4) of Dehnition 



that 



2.3 



It follows 



■ ■ ) ' ii ;; I - '....' 1 ' ' ' " ' .) ■»!!■■ • (4-12) 
Finally, we perform our last application of the fundamental theorem to obtain 

I||<=)'2 +1 (7' t hlMfi _ \\\p l i+ 1 (hi t ^||2 

2 H v l > b ni 2 >\\H x 2" O n ' 9 

= E r / V'ftegjft^ |) : V s 6«+>~,t n) |) fed?. 

By applying the same calculations to the previous equations as those leading to 
d4.11|) , we find that 



- i"r. 



2 



H x 



CAt] e CAt . 



(4.13) 



Combining (4.11|, (4.12), and (4.13) gives 

j;J/ (t,T',W')|| ff A < ||W^( — ,*„,- 



||CJ9+ 1 /+' , /Ml ^ ||o9+ 1 i' ™ j- In MI „2CAt , ^ a+ (JlCl\t , _CAn 



which immediately leads to the conclusion 



N 



\Q^y,u/)\\ H x < ||e^(0,0,0)||^e" CAt + CAtEe eCAt 

0=1 

<||e^(0,0,0)||^e ct '+CT'e ct '. 



(4.14) 



Finally, let us estimate ||0f^(O, 0, 0)|| ff x. By definition, we have that 

At this point, we can apply each of the time-derivatives to $(0, 0,0) and use the 
definition (4.1) to translate time-derivatives into spatial derivatives. The maximal 
number of spatial derivatives for each time derivative is max{a, 1} and consequently 
q + 1 time derivatives translates into at most (q + 1) maxja, 1} spatial derivatives. 
As a consequence, we can conclude that 

lieft 1 ^ ' ' )!!^ = lie^i(o,o,o)|| fffc -(, + ij»«{ a .i} < c|KIU*. 



From this, and (4.14), we conclude that the lemma holds also for \l\ = q+ 1. □ 
Using the previous lemma, we now deduce time regularity of the source term F. 
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Lemma 4.8. Let $ be the Strang splitting solution in the sense of Definition J±.l 
and (4.1). Then, for any 3max{a, 1} < k G N such that \\u$\\ H k < C, we have 

\\v 2 t F(t,T,Lj)\\ Hk - 3inaxia , 1} < c, (t,T,«) e n At . 

Proof. Let I — k — 3max{a, 1}. Let i and j denote any one of t, r, or u>. An 
arbitrary component of V^-F can then be written Fij := didjF. By definition, we 
have that 



F, 



didj 



1 



(tf t + div(i9«(i?))+i? T -A(i?) 



1 



My + drtiij - Afaj) + - div (l?ij«(T?) + ^iV(-dj) + tijVtfi) + 0v(0ij)) ■ 



By applying the Holder inequality and Sobolev embedding, we make the gross 
overestimation (note that we can have £ < 2): 

\\Fij\\H' < |||V?*|| H « + \\V*$\\ He+ « + ||Wy • l>(0)||^ + HtfydivuO?)!!^ 

+ || • ©(fy)^ + H^divu^OH^ + HV^ • v(4i)\\ B * 
+ Wfy divu(^)||^ + ||W • v(ihj)\\H* + ll^divu^ij)!!^ 



\v 2 t n Hl+a + ciiv^n^+x (\\v(#)\\ H z +2 + ii div«(i?)ii^ +2 ) 



+ 2\\V t #\\ fft+ i\\v(Vt#)\\ir*+z +2||Vti?|| H i+2||i»(V t i?)||Hi+ 1 

+ C\\#\\ Hi+3 (\\v(rf0)\\ H i + || div«(V?0)|| H *) 

< C (||V t 3 tf|| H . + ||V^||^ +m a X{Q , 1} + ||V t tf|| H « + • 

Since #(0,0) := uq € Lemma 4.7 tell us that the right-hand side is bounded 
and hence our proof is complete. □ 

The next lemma is the main property giving second order convergence. Observe 
that the result does not depend on our specific choice of operators A and B. 



Lemma 4.9. There holds 



Proof. We will prove Lemma |4.9| by direct calculation. Let us begin by estimating 




FuQf,tn, ). Since F(t,t n , ^f) = 0, we have that 

F (— t 

Similarly, we see that (4.3) yields 



t'n ■ 



..„. J=-[A,Bp). 



F T 



-\A,B]{§). 



(4.15) 



(4.16) 



It only remains to estimate F t (^f,t n , ^f). However, as F(t,t n , ^f) = 0, we must 
have that F t (^f,t n , %-) = 0. This, together with (4.15) and (4.16), concludes the 
proof. □ 

Using the previous lemma, we can now prove that the error produced along the 
diagonal (t/2,t,t/2) is second order in At. 



Lemma 4.10. Let i? be the Strang splitting solution in the sense of Definition J^.l 
and (|4.ip. Then, if u £ H k , 

t 



F 



t 

-,t, - 
2' ' 2 



< C(Aty 



H k-3l 
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Proof. Since f (%-,t n , = 0, a Taylor expansion provides the identity 



F 



t t 
-,t, - 
2' ' 2 




(1-1) (jj< 

i/i;(|)'^(i-i) (|)<i-i> 



ds (4.17) 



where the last equality is an application of Lemma 4.9 By taking the H k 3max{a,i} 



norm on both sides of (4.17) and applying the previous lemma, we gather 

< C(At) 2 , 



.t t, 



which concludes the proof. □ 
Proof of Lemma \4-3\ By performing the same calculations as those in the proof 



of Lemma 3.5 with the new estimate on F given by Lemma 4.10 we obtain the 
estimate 

~d t || e(t)|| H *-3 < C(T) ((At) 2 + \\e(t)\\ Hk - 3n ^ {a , 1} ) , t e (0,T). 

Since e(0) = 0, an application of the Gronwall inequality to the previous inequality 
gives 

||e(t)|| ff *-,»«c-.i} < t (At) 2 C, t e [0,T], 
which concludes the proof of Lemma |4.3| and consequently also Theorem |1.2| □ 



5. Applications 

In the previous sections we have established well-posedness and convergence 



rates for both Godunov and Strang splitting applied to (1.1). In this section we 



will examine a range of different equations that are all of the form ( 1.1 1 with A and 



v being admissible (in the sense of Definition 2.3 1. Our primary goal is to equip the 
reader with some relevant applications of our framework. Let us for the convenience 
of the reader repeat our equation and assumptions here. We are working with the 
equation 

u t +div (uv(u)) = A(u), u\ t =o = u , 



where the operators are required to be admissible in the sense of Definition 2.3 

Given a specific equation, the only non-trivial conditions in Definition |2.3| are 
the well-posedness (6) and the commutator estimate (5). The following lemma is 
of help to determine the latter. 

Lemma 5.1. Let a 6 (0,2) and I be any multi-index. The standard differential 
operator A(u) = D l u and the fractional Laplacian A(u) = (— A) a / 2 u both satisfy 



(3)-(5) in Definition 2.3 



Proof. Conditions (3) and (4) are trivially satisfied. It remains to prove that A 
satisfies the commutator estimate (5). In the case of the fractional Laplacian 
(A(u) = (— A) Q / 2 w), the commutator estimate was proved in Corollary 2.5 in [T7] . 
See also the lecture notes by Constantin [3] . For the standard differential operator, 
the result is immediate from the standard Leibniz rule. □ 
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Remark 5.2. By combining the usual Leibniz rule with the corresponding rule for 
the fractional Laplacian, the previous lemma can be extended to also include oper- 
ators on the form A(u) = (—A) a / 2 D L u, that is, any mix of fractional and standard 
derivatives. 

5.1. Burgers type equations. If we restrict to one spatial dimension (N = 1), 
the only valid velocity operator v satisfying requirement (1) of Definition |2.3| is 



v(u) = au, a ER. 

Thus, the type of equations we can consider consists of a Burgers term and a linear 
differential term. In the literature one can find several equations of this type that 



are well-posed in the sense of (6) in Definition 2.3 Some examples are: 
u t + (u 2 ) x = u xwx , (KdV), 
u t + (u 2 ) x = -(-d 2 ) a l 2 u, (viscous Burgers), 

u t + (u 2 ) x = -u xxx + u xxxxx , (Kawahara). 

For the viscous Burgers equation, a > 1 is required to ensure well-posedness [20 . 

5.2. Quasi-geostrophic flow. The following equation has been proposed as a toy 
model for strongly rotating atmospheric flow 

u t + div(uv(u)) = 0, v(u) = curl(-A)~ 1/2 it, (5.1) 

where u is the potential temperature, v(u) is the fluid velocity and the equation is 
valid in two dimensions (N = 2). The reader can consult [H [7] and the references 



therein for more on the physical aspects of the model. The equation (5.1) is in 
the literature referred to as the quasi-geostrophic equation and has in the recent 
years been the subject of numerous analytical studies. This recent interest was 
probably sparked by the Constantin, Majda, and Tabak, paper [7] in which they 



give numerical evidence for a connection between the blow-up of solutions to (5.1) 
and the three dimensional Euler equations. Though the precise type of blow-up has 
been later dismissed by D. Cordoba [9], there remains hope that understanding the 



behavior of solutions to (5.1 1 can aid in characterizing blow-up of solutions to the 
incompressible Euler equations, a long standing open problem. As a consequence, 
most of the recent studies concerns regularity of solutions to equations of the form 



(5.1). A particularly well-studied case is the dissipative quasi-geostrophic equation 
u t + div{uv(u)) = A(u), v(u) = curl(-A)~ ,3/2 u, (5.2) 
where A(u) = — (— A)"/ 2 u. Through the papers (3] [6j El |T9] , it is has been shown 



that ( |5.2[ ) is well-posed in the sense of (6) in Definition 2.3 for a, (3 G [1, 2]. Since 
in addition divu (it) = 0, (2) in Definition |2.3| is also satisfied. Hence, (5.2) with a, 
P € [1,2] is admissible in our framework. 

5.3. Aggregation equations. The following active scalar equation has been pro- 
posed [23l [24j [25] as a model for the long-range attraction between individuals in 
flocks, schools, or swarms: 

u t + div(uv(u)) = 0. (5.3) 
The velocity is determined as the convolution of u with an interaction potential: 

v(u) — V$ * u. 

In these models, u is the density of individuals and v(u) incorporates the pairwise 
attractive forces between individuals in the flock. Two relevant examples of poten- 
tials $ in applications are the radially symmetric $ = 1 — e' x ' and $ = 1 — e~l x l . 
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In the papers P3 15] , sharp conditions are derived on the potential $ under which 
solutions of (5.3) blows up in finite time. In particular, whenever 



1 



dr < oo, 



solutions of (5.3) collapse to a point at the center of mass in finite time. Thus, the 
relevant case $ = 1 — leads to blow up in finite time. This clearly non-realistic 
behavior can be attributed to the lack of any effect incorporating collision avoidance 
in the model. Since collision avoidance is a short-range phenomena, a simple way 
to incorporate it is to add diffusion to the model 

A(u) = -{-A) a ' 2 u. 



u t + div (u v(u)) — A(u), 



(5.4) 



The equation (5.4) has been the subject of the recent studies [TTJ [T^l H5] and well- 
posedness in the sense of (6) in Definition 2.3 has been established for a € (1, 2]. In 



view of this well-posedness result and Lemma 
provided we can verify condition (2) in Definition 2.3 
By applying the Holder inequality, we see that 



5.1 ( 5.4 ) is included in our framework 



dxvv{f)4> dx 



A<f>(y)f(x-y)<p(x) dxdy < \\f\\ L 4<P\\ 



LP 



|A$| 



Hence, div v(f) satisfies (2) in Definition 2.3 for any potential $ satisfying A<f> e L 1 



Comparing this condition with the results in the paper |21j . we see that we can 
include all the potentials for which well-posedness is known. For instance, the case 
$ = I — e~l x l , since 



A$ = div 



e -|x| =_ e -N 



N- I 



-M 



x^O, 



and 



N -1 



dx <C, N>2. 



{M<i} 



5.4. Magneto geostrophic dynamics. In the recent paper [2], well-posedness 
was established for the following class of active scalar equations 



Ut + div(uv(u)) — A(u), v(u) = divT(u), 
where T is a matrix of Calderon-Zygmund operators satisfying 

divu(u) = divdivT(u) = 0, 



(5.5) 



and A is the Laplace operator A(u) = Au. The equation (5.5) can be seen as a 
generalization of the quasi-geostrophic equation (5.2) (with a = 2) since the latter 



can be obtained from (5.5) with a particular choice of T in 2D. The equation (5.5) 



is also included in the framework considered by Constantin in [5] . 

The physical motivation for the study [13] was that (5.5) with a particular choice 
of T has been proposed as a model for magnetostrophic turbulence in the Earth's 
fluid core. See the paper [141 and the references therein for more on this application. 
Since (5.5) is well-posed in the sense of (6) in Definition 2.3 our framework does 
indeed include this class of active scalar equations. 



Appendix A. Proof of Lemmas 12.11 and 12.2 



In this appendix we have gathered the proofs of Lemmas 2.1 and |2.2[ Both 
lemmas have been used in an essential fashion throughout the convergence analysis. 
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Lemma 2.1. Let k > 6 and v be an operator that satisfies Definition \2.3\ Then, 



£ / V fl (div(/«(/))) : V s / dx 



<C\\f\\„*-4f\\H», f^H k . 



s=0 



Proof. Let us confine to the three dimensional case (N = 3) as the other cases are 
almost identical. The product rule provides us with the identity 



T [ V s dw(fv(f)) : V'fdx 

V( / V s (V/ • «(/)) : V s f dx + [ V s (fdWv(f)):V s fdx 

k 

= £|/fl + l3l- 



< 



In the remaining parts of the proof, our strategy is to bound each of If and Z|, 
s = 0, . . . , k, separately. We begin with Ij 5 . 

1. By applying the Leibniz rule to 1^ (with multi-index notation a = (o^ , a2, 0:3)), 
we obtain the following expression 



P i= E / V a (Vf-v(f))V a fdx 



|o|=s 



Ql Q2 Q3 



££££ 

|o|=» ii=0 i2=0 23=0 



Oil \ I Oti \ ( CK3 

*1 



V „ ; - .. - - ■ V 



«2 / \ *3 

Qs-i!-i 2 -i3 f 



dx^dy^dz^ \dx ai ~ il dy a2 ^ i2 dz a ^'^ J J dx^dy^dz 013 



d s f 



(2-1) 



dx. 



Let us now consider four separate cases of i\ + ii + 13 in the above quadruple sum. 

(i) If i\ + %2 + «3 = k = s, i.e., (ai, a2, 03) = (ii, 12, *3)> the above term can be 
rewritten as follows 



gkf 



«(/) 



dx ai dy a2 dz a 3 J dx ai dy a2 dz a ^ 



d k f 



1. 



P: .v 



d k f 



dx a ^dy a2 dz a 3 



v(f) dx 



dx 
1 

v 2 



gkf 



divv(f) dx 



< 



d k f 



dx ai dy a2 dz c 



dx a ^dy a2 dz a 3 

diV«(/)||icc <C||/||^||«(/)||K3<C||/||^||/|| Hfc - 2 



L 2 



where we have used the Sobolev embedding H 2 C L°° and fc > 6. 
(ii) If 2 < + «2 + *3 < ft — 3, we find 



QS- 11—12—13 J 



R « V dx^dy^dz* 3 \dx a ^-^dy a ^-^dz a3 - i3 J J dx ai dy a2 dz a3 



dx 



< 



Qil+i-2+i3 f 



dx^dy %2 dz^ 



d 



s— 21—22—23 



/ 



Q x a 1 -iiQya 2 -i2Q z ci3-i3 



L' 1 



d s f 



dx ai dy a2 dz a3 



<tf||/IU*l|/IU*-»ll/IU»- 
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(iii) If k — 2 < ii + i% + is < k — 1, we find 

Qs-ij_-i 2 — i 3 j 



< 



R « \ dx^dy^dz^ V 

Qil+i2+is. j 



d x ai-iiQya 2 -i2d z a 3 -i 3 J J dx ai 8y a2 dz 



dx 



dx il dy l2 dz i3 



L 2 







S — 1\— 12— 13 



f 



dx a i~ l idy a2 ~ l2 dz a 3- 



d s f 



dx ai dy a2 dz a 3 



L 2 



<C\\f\\ H 4f\\ H 4f\\ Hh <C||/||h*-2||/||^. 



To conclude the second last inequality, we have used that k — 2 < s < k and hence 
that s — i\ — %2 — «3 < 2. The inequality then follows from the embedding H 2 C L°° . 
(iv) If < i\ + 12 + «3 < 1, we find 



< 



Qii+'h+ia j 

r n V. " dx ll dy i2 dz is 

Qil+l2+i 3 j 

dx^dy %2 dz i3 







S — H— 12 — «3 



Q^-HQyOii-izOzaa-is 
qs— ii— i 2 — is j 

dx ai - ll dy a2 - l2 dz a3 - i3 



dx 



dx a ^dy a2 dz a 3 

d s f 
dx ai dy a2 dz a3 



L 2 



L 2 



<c\\f\\ H 4ft H * <c||/|| fl *- a ||/| 



Hence, by applying (i)-(iv) in (2.1|, we see that 

k 

52\Il\<C\\f\\ Hk -4f\\ 2 Hh . (2.2) 

s=0 

2. We now bound the I| terms. The standard Leibniz rule provides the identity 
I '=Y1 I V a (fdwv(f))V a f dx 

(2.3) 



l«l= 



Oil Oi2 «3 



EEEE 

|a|=s ii=0 12=0 i 3 =0 

/ Qil+i2+is j 



oti\ ( a- 2 \ ( Oi-i 
h 



R « \dx ll dy l2 dz i 



div v 



h ) V h 



Qxax-hQyaz^Qzas-is J J dx ai dy a2 dz° 3 



dx. 



Let us again consider four separate cases of i\ + i 2 + 23. 

(i) If i\ + i 2 + 13 = 0, we apply the Holder inequality and (3) in Definition 2.3 
which yields to deduce 



/divw 



d s f 



< 11/11 



div v 



dx ai dy a2 dz a 3 J J dx ai dy a2 dz a3 
d s f 



dx 



dx a ^dy a2 dz a3 
(ii) If 3 < ii + 12 + «3 < k — 2, we find 



L 2 



ii/Hh. <c\\f\\ H >-4fr H >. 



< 



Qil+i2+i3 j 

dx^dy i2 dz l 3 



diw 



QS— il— i 2 — 13 f 



Q x ai-hQya 2 -i 2 g z a 3 -i 3 J J dx ai dy° 2 dz a 



d s f 



dx 



dx^dy^dz* 3 



div v 



Qs—ix—i2—i 3 f 



Q x ai-h Qya 2 -i2d z a 3 -i 3 



d s f 



dx ai dy a2 dz a 3 



< C\\f\\ Hk -2\\v(f)\\ H 4f\\ HS < C||/|| H »_ a ||/||^. 
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(iii) If 1 < i\ + i-2 + 13 < 2, we find 



< 



dx^dy^dz^ 

Qil+i2+i:i j 



divt) 



Qs— ii— 12— is, j 



Q x a 1 -i 1 Qya2-i2Q z a :i -i: i J J dx ai dy° 2 dz 



d s f 



dx 



d x ilQyi2Q z i3 



div v 



Qs-i!-i 2 -i3 j 



dx^-^dy^-^dz^"^ 



L 2 



d s f 



dx ai dy a2 dz a 3 



< C||/||H«»||t;(/)|| ff .||/|| ff - < C\\f\\ H >-4f\\ 2 Hh . 
(iv) If k — 1 < i\ + 12 + «3 < k, we find 



< 



^1+12+13 f 

dx^dy^dz^ 

Qii+i2+'ki j 



divt) 



Qs-i 1 -i 2 -i3 f 



dx^dy^dz^ 



L 2 



div v 



Q x a 1 -i 1 Qya 2 -i2Q z a3-i3 

QS — 11—12 — 43 f 



d s f 



dx OL1 - il dy ot2 - i2 dz a z- i3 



dx ai dy a2 dz a z 
d s f 



dx 



dx ai dy a2 dz c 



< C\\f\\ Hk \\v(f)\\ H 4f\\ H s < C\\f\\ Hk - 2 \\f\\ 2 Hk 
Applying (i)-(iv) in (|2.3|) gives 



k 

E 

s=0 



m<c\\f\\ Hk -4f\\ 2 Hk . 



Together with (2.2), this brings our proof to an end. 



Lemma 2.2. Let k > 4. Then the following estimates hold 

<C\\g\\ H 4ft H *, 



k 

E 

s=0 



j V s div (/«(<?)) : V s / dx 
£ [ V s dW(gv(f)):V s f dx 



<C|| 5 ||^ + l||/|| 2 ^. 



□ 



(2.4) 



(2.5) 



Proof. The proof of (2.4) is easily obtained by the calculations of the previous 
proof. To prove (2.5), it is only step (i) in part 1 of the previous proof that does 

aake t 

v(f) 



not go through. Instead, we now make the calculation 

dx ai dy a2 dz a 3 



dx ai dy a2 dz a 

d k g 

dx ai dy a2 dz a 3 

and the proof is complete. 



dx 



< 



Hf)W L ~ \\f\\ H « <c\\f \\W\\ f\\ 



H 2 i 



L 2 



□ 
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